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Changes in Dark Matter Properties After Freeze-Out 
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The properties of the dark matter that determine its thermal relic abundance can be very different 
from the dark matter properties today. We investigate this possibility by coupling a dark matter 
sector to a scalar that undergoes a phase transition after the dark matter freezes out. If the value 
of Adm h 2 calculated from parameters measured at colliders and by direct and indirect detection 
experiments does not match the astrophysically observed value, a novel cosmology of this type could 
provide the explanation. This mechanism also has the potential to account for the "boost factor" 
required to explain the PAMELA data. 

PACS numbers: 95.35.+d 



I. 



INTRODUCTION 



The amount of dark matter (DM) in the universe has 
been measured precisely by astrophysical and cosmolog- 
ical probes: (0 DM h 2 ) astro = 0.106 ±0.008 Q. The lead- 
ing candidate for this DM is a new weak-scale neutral 
particle. If the universe follows a standard thermal his- 
tory, the DM density can be derived from measurements 
of the properties of this particle. The crucial input is the 
particle's thermally averaged annihilation cross section, 
(a a v). When the annihilation rate becomes too slow to 
keep pace with the expansion of the universe the DM 
particle freezes out, leaving behind a relic density of DM 
that merely dilutes as the universe expands. 

With the turn-on of the Large Hadron Collider (LHC) 
and a host of direct and indirect detection experiments 
coming on-line, there is hope that the nature of the 
DM particle will be measured thoroughly enough that 
(cr a t>) can be computed. Then a prediction of the ther- 
mal relic abundance, (Odm h 2 ) part i C ie, can be made. If 
(0 DM h 2 ) part icie = (^dm h 2 ) astro , this will be strong ev- 
idence that the universe has a standard thermal his- 
tory back to the DM freeze-out temperature, Tf a (typ- 
ically tens of GeV for weak-scale DM). This would ex- 
tend the successful predictions of Big Bang Nucleosyn- 
thesis (BBN), which demonstrate a thermal history of 
the universe only back to temperatures of several MeV. 

On the other hand, if the calculated relic density 
does not equal the measured one, this will be ev- 
idence for physics beyond minimal thermal DM. If 
(0 D m h 2 ) par ti C i e < (Odm h 2 ) astro , it is possible that we 
have not identified the dominant source of DM or the 
DM was produced non-thermally as a decay product of 
another particle 0. Conversely, if (Odm h 2 ) par u c ie > 
(Odm h 2 ) astro the thermal relic abundance of the DM 
must have been diluted, perhaps by a late production of 
entropy [|[ or a modification of the expansion history of 
the universe [3, In the present work, we explore a 
novel possibility that can obtain either direction of this 
inequality: a change in the properties of the DM itself 
between Tf a and the present. Time-dependent DM has 
been considered in another context in attempts to relate 
DM and Dark Energy @, @] . 



Relevant changes in the attributes of the DM particle 
can occur if there is a field whose vacuum expectation 
value (VEV) changes during the crucial epoch between 
Tf and BBN. If this field influences the mass or cou- 
plings of the DM particle, there can be a dramatic effect 
on the relic abundance one would calculate based on the 
properties of the DM particle measured today. Here we 
present a simple model that illustrates how this mecha- 
nism could be realized. We discuss some constraints on 
scenarios of this type, and we study the phenomenology 
that should accompany the late-time phase transitions 
typical of this class of models. 



II. 



A LATE-TIME PHASE TRANSITION 



To change the DM properties, we suppose there is a 
phase transition (PT) after Tf Q 8]. In the model con- 
sidered here, this PT occurs in a new sector containing 
a Standard Model (SM) singlet P. We couple the PT 
sector to a model for the DM in the next section. The 
PT will modify both the mass and couplings of the DM 
particle in this model. 

Rather than introducing a new field P, one might 
instead try to modify the properties of the DM after 
freeze-out via the electroweak PT. This does not work 
for electroweak-mass DM using the minimal SM Higgs 
phase transition Q . Unless the dynamics of this PT are 
modified (or the initial DM mass is very large), the tem- 
perature of the PT is typically greater than Tf a , and 
the DM properties would not be modified between Tj a 
and the present day. On the other hand, if the Higgs 
boson sector is non-minimal, it is possible that the elec- 
troweak transition temperature might be lowered sub- 
stantially (see e.g. [Tot]). 

The new singlet field P is initially stabilized at the ori- 
gin in the early universe by a thermal mass term [Til . [T3 | . 
As the universe cools, P undergoes a PT at a temper- 
ature Tpx < Tf a ~ todm/20, and develops a non-zero 
VEV, (P) = vp. For the PT to have a significant effect 
on the DM properties (perhaps by generating a large ex- 
cursion in the DM mass A m = \bm-p vp) typically re- 
quires Vp 3> TpT [H|. 
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We take the potential for P to be 

V P (T = 0) = -±\m P \ 2 P 2 + ±P\ (1) 
which induces 



\mp\ 


A 


vp 






TpT 




D 


4.0 GeV 


1.5 xicr 5 


2.5 TeV 


0.10 


9 


33 GeV 


40 GeV 


0.77 



TABLE I: Benchmark phase transition parameters. 



v P (T = 0) = V6|mp| 2 M 



(2) 



below Tpt- The Z2 symmetry of this potential (P — > 
— P) means there is a danger of forming domain walls. 
We can retain the form of the potential while avoiding 
domain walls by softly breaking the I2 with a very small 
cubic term, making this symmetry only approximate [1 31 ] . 

A large hierarchy between vp and Tpx in this scenario 
requires that the coupling responsible for inducing a ther- 
mal mass for P be considerably larger than A. This can 
arise if P couples to other states that are approximately 
massless when vp = 0. Such states can emerge if P is 
part of a larger "hidden" sect or, p erhaps coupled to the 
SM only via a "Higgs portal" For concreteness, 

we consider additional fcrmionic fields coupling to P ac- 
cording to C 3 XpQiPQi Qi- Since the Q's have no other 
mass terms (which would violate the Z2 of P) , these cou- 
plings contribute to the temperature-dependent mass of 
the P field, strongly trapping it at the origin. When vp 
shifts to its non-zero value, the Q's acquire a mass of 
ApQ; Vp, typically of order a few hundred GeV. 

At high temperatures and near the origin of P, the 
potential is approximately [HI, [l2j 

Vp{T) = -\{\m P \ 2 - ^\% Q T*)P 2 + iAP 4 , (3) 

where Xpq is the (universal) coupling between P and 
the Q's and Nq is the number of Dirac Q fields. This 
potential gives a PT temperature of 



Tpt 



NqX 2 



(4) 



PQ 



ally 
17. 



Strong trapping of the P field at the origin typical! 
leads to a brief period of thermal inflation (TI) [1 
[l8j . The vacuum energy density during TI is p vac = 
\mp\ 2 vp/4:. If TI ends at Tpt by the instantaneous decay 
of the P field to radiation, we can estimate the reheating 
temperature Tpjj via conservation of energy: 



rrii 

1 RH 



45|m P | 4 
,g? H it 2 A 



36 g 



PT 1 |4 
TOP 



.9 



RH 



Nq 2 \ pq : 



(5) 



where c/» is the effective number of relativistic degrees 
of freedom. Reheating can dilute the DM abundance. 
Although this is not the dominant effect that we wish 
to explore, it can be of quantitative importance. This is 
also the reason why we rely on thermal corrections, rather 
than an additional cubic term in the tree-level potential, 
to trap P at the origin. With a cubic term, the trapping 
need not turn off as the universe supercools and could 
lead to a severe dilution of the DM abundance. 



To estimate this dilution, we first assume there are no 
new sources of entropy during TI. This fixes rif /sf — 
ripx/spx, where n and s are the number density of the 
DM and entropy density of the universe respectively. No 
DM is produced in the reheating process, implying n^T = 
«rh- Once TI ends and reheating completes, the new 
conserved quantity is tirh/srh = "pt / srh- The dilution 
factor, D, is 



"PT _ SPT Ufp 

srh srh Sf 



I* J PT \ n fo 



RH T 3 



RH/ S fo 



■UJL = D y,"ltL. ( 6 ) 



Sfo 



Taking into account the change in the mass of the parti- 
cle, the present abundance is given by 



(f^DM h 2 ) ast ro = D X 



'dm 

v P =0 
DM . 



qup=0.2 
DM ' 



(7) 



This can differ dramatically from (Hdm h 2 ) part i C ie, as we 
will see in the next section. 

In Table [J we exhibit a benchmark point that gives a 
first-order PT with a transition temperature Tpx <C vp. 
To obtain this feature, the value of A is small. This in- 
teraction obtains additive corrections of the form AA = 
5^(cb A 2 — c/ Aj)/(167r 2 ), where the sum runs over bosons 
and fermions that couple to P, and the Cj are 0(1) coef- 
ficients. For the benchmark couplings, the small value of 
A is technically natural. 

The value of Tpx for the benchmark point is also large 
and could exceed a typical value of Tf Q unless the mass 
of the DM particle is many hundreds of GeV. Smaller 
values of Tpt can be achieved by reducing the value of 
|mp| 2 . This leads to light excitations of P that can be 
phenomenologically problematic - it is difficult to make 
them decay quickly enough to avoid BBN constraints 
while not disturbing the evolution of supernovae. 



III. 



A DARK MATTER SECTOR 



There are many possibilities for the DM sector, all of 
which could work with the generic phase transition mod- 
ule we presented in the previous section. The particular 
DM sector we consider is a "level-changing" model, con- 
sisting of three fermions with the same quantum numbers 
as the Higgsinos and Bino of the minimal supersymmet- 
ric SM: a vector-like pair of SU(2)l doublets tpL and ipi 
with the appropriate hypercharges, and a gauge singlet 
ips ■ All fields in this DM sector are charged under an ex- 
act X — > —X symmetry (independent of the approximate 
Z2 of P), implying that the lightest of these particles is 
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I 1 


Ms 


A s 


Ai 


A 2 


1.3 TeV 


0.68 TeV 


-0.070 


0.020 


0.010 


niDM^p = 0) 


itidm(«p 7^ 0) 


T fo (v P = 0) 


T f o(vp + 0) 


1.3 TeV 


1.0 TeV 


65 GeV 


52 GeV 



TABLE 

(S!dm h 2 ) pa 



II: 

rticle 



Benchmark 



parameters realizing 



absolutely stable. The DM sector Lagrangian is 

C 3 nipL-ip L + \ 1 H-ip L ip s +\2H*-ip L <ip s (8) 
+ (fj, s + X s P) ip s ip s + h.c, 

where H=(G+, -±(H° + i G°)) T is the SM Higgs boson. 
The resulting "neutralino" mass matrix is 



M° = 





/'- 







(9) 



with electroweak VEV (H°) = v„ = 246 GeV [H. 

Within this model, it is not difficult to obtain 
(Qdm h 2 ) P articie > (^dm h?)astro- As an example, we 
consider the benchmark parameter point given in Ta- 
bles U and [Til At high temperatures vp = 0. There 
the DM is a nearly pure combination of the doublets 
ip L and Vl : x ° ~ l/V2ipL + l/y/2ipz + eip s , with 
e ~ (Ai - A2)wff/(4^ s — 2^). The thermal relic abun- 
dance of this state is nearly identical to that of a pure 
Higgsino. This is set by its annihilation to pairs of W 
bosons, and is given by [l^| 



n Vp=0 h 2 
"dm '* 



0.1 



todm 



ITeV/ 



V 



(10) 



including coannihilation with the heavier "charginos" . 

The mass and composition of the DM change after the 
PT. For the parameters in the Tables, the lightest of the 
DM-sector particles is nearly pure singlet post-PT. Using 
Eq. ([7]), its relic density is (Odm h 2 ) as tro = 0.1. This is 
the value measured by astrophysical probes. However, it 
is considerably different from the value one would recon- 
struct from measurements of the DM particle Lagrangian 
today, assuming one measured the relevant couplings but 
did not take into account the non-canonical cosmological 
effect described here [33 1. 

The dominant contribution to the apparent particle an- 
nihilation cross section, assuming the relevant particles 
and their couplings can be measured, is the s-channel 
exchange of a P going into QQ. Assuming a standard 
thermal history, the predicted relic density is approxi- 
mately given by 



(f^DM h ) pa rticle 
yielding (17 DM h 2 ) p 



0.02 



n q (x P q x s y 



G 



m D M 



V 



TeV/ 



(11) 



Jparticie — 45 for the benchmark, more 
then two orders of magnitude larger than (JIdm h 2 ) astr o- 



Even if the PT-sector particles are not discovered at col- 
liders, the properties of the DM today will differ from 
those at freeze-out. These properties can potentially still 
be deduced by direct and indirect detection searches for 
DM. 

We obtained (Odm h 2 ) partlc i e > (Qdm h 2 ) astro in this 
example. A different choice of mass matrix (Eq. ([9|)) 
can lead to the opposite relationship. When this is the 
case, the value of (a a v) should increase after the PT, and 
the DM can potentially recouple after thermal inflation. 
Demanding that the DM stay frozen out gives a bound on 
the allowed change in the relic density. Non-recoupling 
of the DM after reheating requires 



np T (a a v) v ^° < 1.66(<£ 



rp2 

RH^l/2 J R.H 

Mpi • 



(12) 



where Mpi is the Planck mass. A similar condition holds 
for the initial (vp — 0) freeze-out cross section and tem- 
perature. Combining these expressions and accounting 
for redshift from freeze-out to the PT gives 



(o a v 



v P =0 



< 



<j 



PT 



T 2 T 



v P =0 - 



(13) 



Here gl P is the effective number of relativistic degrees 
of freedom calculated at TJ P=0 . Using the standard ap- 
proximate solution to the Boltzmann equation [20] to re- 
late (<7 a v) to SIdm/i 2 , along with Eq. (0, leads to the 
constraint 



(Odm h 2 ) pa rticle 
(Odm h 2 ) 

astro 



> 



fJ 



RH 



?RH 



(14) 



1 fo 



A large change in the apparent relic density without 

^ j 

fo 



recoupling requires a hierarchy between Trh and T^^ (l 



To avoid disturbing BBN, Trh must be larger than about 
10 MeV. Taking a typical Tf a of tens of GeV, the appar- 
ent relic density can be reduced by a factor of a thousand. 
In practice we find it difficult to obtain such low reheat- 
ing temperatures simultaneous with the large vp needed 
to make a significant shift in the DM properties. 



IV. PHENOMENOLOGY OF A LATE PHASE 
TRANSITION 

For the PT to happen after DM freeze-out, the mass 
of the physical P excitation ~ \rnp\ should be light. The 
existence of a light P is the most generic feature of the 
mechanism presented here, and so it is worth considering 
its phenomenology in some detail. The symmetries of the 
model allow the Lagrangian term C 3 (Xph /2)P 2 \H\ 2 , 
coupling P with the SM Higgs boson. The resultant mix- 
ing with the Higgs boson gives two mass eigenstates, p° 
and h . The mixing angle is given by 



tan 26 



6 Xp H vp vh 
X H v 2 H ~ Xv 2 p ' 



(15) 



4 



where C 3 X H (H°) 4 /4\ 

The relevant phenomenological constraints and signals 
depend on the precise mass of the p°, which in principle 
could range from tens of GeV all the way down to a 
fraction of an MeV. Mixing allows the p° to be produced 
in association with a Z°, or to appear in meson decays. 
For TOpO < 100 MeV, astrophysical constraints similar to 
those for axions (2l| become important. 

For DM masses near the weak scale, the natural value 
of the p° mass is on the order of a few GeV. For the 
parameters in Table U and a moderate mixing angle, 
TOpO ps 6 GeV. In this mass range, the p° could be pro- 
duced in Upsilon (T) decays. To lowest order (22| . 



r(T^p° 7 ) 



sin 



r(T 



) 



l-^f • (16) 



Requiring BR(T -> p° 7) x BR(jj° 



T T 



') < 10" 



gives a modest bound on the mixing angle of 9 < 0.3 for 
BR(p° — * t + r~) = 1. For a 6 GeV p°, decays to charm 
quarks actually exceed those to r's by a factor of 2 (un- 
less a more complicated Higgs sector allows for a tan j3 en- 
hanced p° couplings to down- type fermions). In this mass 
range, a comparable bound exists from non-observation 
of Z a p°, which would have been seen in Z° h° searches 
at LEP [2J|. For higher masses, top > 10 GeV, the 
bound on the mixing angle strengthens due to the LEP 
constraint: 9 < 0.14. The p° decays are very prompt; 
the lifetime of p° is 2 x 10~ 19 sec for 9 = 0.14. The p° 
branching ratios are identical to a SM Higgs boson of the 
same mass. If the p a mass falls below the B-meson mass, 
bounds on the mixing angle from b — > s P processes [l[ 
are strong: 9~^ 10~ 4 . 

These considerations also provide a way to observe the 
p° state at colliders. For lighter masses (m p o < 8 GeV) 
and large mixing, searches for the rare decay T — > 7p° — > 
7t + t~ may be useful. If the h° is not so heavy that it 
decays to W bosons, then h° — ► p° p° need only compete 
with h° — ► bb [25| . The ratio of the widths is given by 



r( h ° -> gV) = 3 {V^vpVHf Ko - 4mg ) 1 / 2 
T{h°^bb) m 2 (m 2 h0 -4mf) 3 / 2 

with the effective coupling 



to 



Ah v h 2 

■ to 



A vp 
he A g v% 
18 ^ Xv% 



(17) 



(18) 



1 - 2fc 



V H 



where tg = tan# and eg = cos 9. For 9 saturating 
the LEP bound, BR(h° — > p°p°) can approach 40% for 



v P 



> 750 GeV. Then h° Z° 



.0 ^0 70 



p'p 



ibZ° might 



be observable at the LHC if 6-tagging efficiencies are suf- 
ficiently high 26], though it will be challenging. 

Thus far we have not mentioned the decay of the Q's. 
This can proceed via higher-dimension operators. Alter- 
nately, the Q's can decay to quarks through renormaliz- 
able operators if they are SU (3) c triplets and are allowed 



a very small mixing with the quarks of the SM. For the 
benchmark parameters in Table HI vtiq = 250 GeV. Then 
given the latter scenario there is the possibility of pro- 
ducing the Q's directly at the LHC. 



V. DISCUSSION 

The DM model presented here represents an existence 
proof of a general mechanism: DM properties can change 
after freeze-out. We have focused our attention on situ- 
ations where a shifting VEV causes a change in the DM 
mass and composition. A similar effect could occur if the 
coupling that sets the relic abundance of the DM is a 
function of a light modulus. 

In another example, the DM mass might shift so that 
2 todm is approximately resonant with some other state 
in the theory, such as a Higgs boson. With the DM now 
sitting on resonance, one would calculate a tiny thermal 
relic abundance. To implement this scenario using a SM 
Higgs boson is difficult. For a "natural" PT, ttt-dm ~ 
TeV. To access the Higgs resonance, m^o ~ 2 todm ~ 
O(TeV). This implies T h o will be too large to generate a 
strong resonant enhancement of (a a v) . In the presence of 
heavy but narrow resonances, this is a viable mechanism. 

Alternately, the DM itself could remain unchanged, 
but the properties of particles crucial for setting the ther- 
mal relic abundance are modified by the cosmology. Con- 
sider a coannihilating particle, C, nearly degenerate with 
the DM. If the mass of C shifts between Tf Q and now, 
the importance of coannihilation would not be evident 
from low-temperature measurements, and the calculated 
(f^DM h 2 ) P articie would differ from the true value. 

If (FIdm h 2 ) particle ^ (^dm h 2 )astro, it is possible that 
the relic abundance of the DM is actually thermal, but an 
alternate cosmology has altered the DM properties since 
freeze-out. These scenarios are naturally realized if there 
is a light modulus that undergoes a late PT. The field 
responsible for the late time transition, P, could show up 
in future experiments. One possibility is via Higgs boson 
decays: h° — > p° p°. If p° is light enough, it could also be 
produced in rare meson decays. Embedding a model of 
this type in an extension of the minimal supersymmetric 
SM is a direction for future investigation. 

Recent preliminary data from the PAMELA [27l ]. 
ATIC 28], and PPB-BETS [H| experiments report sig- 
nificant excesses of cosmic ray positrons and electrons 
above the expected astrophysical background. This ex- 
cess could be the result of dark matter annihilation in 
our galaxy. However, such a dark matter interpretation 
of these results requires a DM annihilation cross-section 
well above the value that would generate the observed 
dark matter relic density [30]. The cosmology discussed 
here offers the possibility of explaining these indirect 
DM signals while maintaining a fairly standard thermal 
freeze-out picture. What is needed is a DM annihilation 
cross-section that increases significantly between freeze- 
out and today. 
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